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ABSTRACT 



Using a high resolution hydrodynamical cosmological simulation of the for- 
mation of a massive spheroidal galaxy we show that elliptical galaxies can be very 
compact and massive at high redshift in agreement with recent observations. Ac- 
cretion of stripped in-falling stellar material increases the size of the system with 
time and the central concentration is reduced by dynamical friction of the surviv- 
ing stellar cores. In a specific case of a spheroidal galaxy with a final stellar mass 
of 1.5 x 1O 11 M we find that the effective radius r e increases from 0.7±0.2 kpc at 
z = 3 to r e = 2.4 ± 0.4 kpc at z = with a concomitant decrease in the effective 
density of an order of magnitude and a decrease of the central velocity dispersion 
by approximately 20% over this time interval. A simple argument based on the 
virial theorem shows that during the accretion of weakly bound material (minor 
mergers) the radius can increase as the square of the mass in contrast to the 
usual linear rate of increase for major mergers. By undergoing minor mergers 
compact high redshift spheroids can evolve into present-day systems with sizes 
and concentrations similar to observed local ellipticals. This indicates that minor 
mergers may be the main driver for the late evolution of sizes and densities of 
early-type galaxies. 

Subject headings: galaxies: elliptical - galaxies: interaction- galaxies: structure 
- galaxies: evolution - methods: numerical 
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Introduction 



There is recent observational evidence that a significant fraction of massive evolved 
spheroidal stellar systems is already in place at redshift z > • 2. H owever, only a small percent- 
age of these galaxies is fully assembled (Ivan Dokkum et al.ll2008l ) . The galaxies a re smaller by 



a fac t or of three to five compared to present - day ellipticals at similar masses (IDaddi et al. 



20051 : lLonghetti et al.l 120071 : 



Toft et all 120071 : iTruiillo et all 120071 : Ivan Dokkum et al 



2008 
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Cimatti et al. 2008; Saracco et a 



one order of magnitude higher ( 



2008). Their effective stellar mass densities are at least 



van Dokkum et al.l 120081 ; iDamjanov et al.l 120091) with sig- 



nifica ntly higher surface brightnesses compared to their low redshift analogs (jCimatti et al. 
2008h . 



These observations are difficult to reconcile with some current idealized formation sce- 
narios for elliptical galaxies. A simple conclusion from the data is that most early-type 
galaxies can neither have fully formed in a simple monolithic collapse nor a binary merger 
of gas-rich disks at high redshift, unless their increase in size can be explained by secular 
processes such as adiabatic expansion driven by stellar mass loss and/or strong feedback 
( IDamjanov et al.ll2009l ; iFan et al.ll2008l ). Additionally, simple pass ive evolution of the stella r 
population is in contradiction with observations of local ellipticals (Ivan Dokkum et al.ll2008l ). 



Dry (i.e. gas-poor, collisionless) mergers and stellar accretion events ar e the prime 



cand i dates for the strong mass and size evolu ti on of stellar spheroids at z < 2 (INaab et al. 



20061; IKhochfar & Silklbood: bell et dJ lioOfoK iRuszkowski & Springellliooi Iffopkins et al 



20091 ; Ivan der Wei et al.l 120091 ) as the additional presence of a dissipative c omponent in a 



major merger event would limit the size increase (see e.g. ICiotti et al.l 120071 ) . The observed 
ellipticals are already very massive at high redshift, thus we expect from the shape of the 
mass function that minor mergers should be much more common than major mergers until 
z=0. 

Massive early- type galaxies m ay undergo not more than one major merger ( with typi 
cally l ow cool g as content, see also IWhitaker fc van Dokkumll2008l ) since z = 0.7 (IBell et al. 



2006al . see also iMcIntosh et al.l 120081) wi th a significant contribution from minor mergers 
for the mass buildup ( Bundy et al.l 2009 ). The low number of observed major early-type 
mergers is also supported by theoretical evidence that massive (~ 10 n ~ 12 M Q ) halos at z = 2 
typically experi ence only one major merger or less until z = and minor mergers are much 
more common (IGenel et al.l l2008bl ; IKhochfar fc Silk! 120081 ) . On average, this is not enough 
to account for the required mass and size growth (see also Damjanov et ali 2009 ) as ma- 
jor dry mergers at most increase the size of a simple one component sy stem by a factor of 



two and allowing for dark matter halos re duces the size growth further (INipoti et al.l 12003 
Naab & TruuUolbood : Iflopkins et al~1l2008h . 



In this Letter we use, as a proof of principle, a very high resolution cosmological simula- 
tion of the formation of a spheroid with no major mergers below z = 3 in combination with 
simple scaling relations to show that the observed rapid size growth and density evolution 
of spheroidal galaxies can be explained by minor mergers and small accretion events. The 
problem is computationally very expensive. At high redshift the observed ellipticals have 
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half-mass sizes of < 1 kpc with accretin g subsystems of even s maller size. As we know 
from isolated merger simulations (see e.g. iNaab fc Burkertl 120031 ) . to resolve such a system 
reasonably well we require a force softening of 10% of the effective radius, which in our case 
is of the order of 100 pc and we require particle numbers of ~ 10 7 to simulate the galaxy in 
a full cosmological context over a Hubble time. Finally, to accurately follow the kinematics 
high force and integration accuracy are required. 



2. Size evolution and minor mergers 



Using the virial theorem we make a simple estimate of how an initial one-component 
stellar systems evolves when mass in stellar systems is added. We assume that a compact 
initial stellar system has formed dissipatively from stars. This system has a total energy 
Ei, a mass M», a gravitational r adius r 9ii , and the mean square speed of the stars is (vf). 
According to the virial theorem (IBinney fc Tremaind 120081 ) the total energy of the system is 



Ei = K i + W i = -K i = -W i 

1 , . 2 , 1 GMf 
= — Mi(vf) = -. 

2 2 r 9ii 



(i; 



We then assume that systems are accreted with energies totaling E a , masses totaling 
M a , gravitational radii r a i and mean square speeds averaging (v%). We define the fractional 
mass increase from all the accreted material rj = M a /Mi and the total kinetic energy of the 
material as K a = (l/2)M a {v^), further defining e = (vl)l(vj). Assuming energ y conservation 
(orbital parameters from cosmological simulations (IKhochfar fc Burkertl 120061 ) indicate that 
most halos merge on parabolic orbits), the total energy of the final system is 



E f = Ei + E a = ~Mi{v1) - ^M a (vl) 
= - l -M l {v}){l + tr 1 ) 

= -\M f (vj). (2) 

The mass of the final system is Mf = Mi + M a = (1 + rf)Mi. Therefore the ratio of the 
final to initial mean square speeds is 
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(vj) _ (1 + rye) 

Similarly, the ratio of the final to initial gravitational radius is 



(3) 



r gJ 



g,i 



[l + rie) 



(4) 



and for the ratio of the densities we get 



Pi 



(l + r/) 5 



(5) 



If during one or more mergers the initial stellar system increases its mass by a factor 
of two then r\ — 1. This mass increase can be caused by one equal- mass merger in which 
case the mean square velocities of the two systems are identical and remain unchanged in 
the final system (Eqn. [3]). The radius incr eases by a factor of t wo (Eqn. HJ) and the density 
drops by a factor of four (Eqn. [5]) (see also lHopkins et al.ll2009[ ). If, however, the total mass 
increase by a factor of two is caused by accretion of very small systems with (v%) « (vf) 
or e << 1, then the mean square velocities are reduced by a factor two, the radius is four 
times lar ger and the den s ity is reduced by a factor of 32 with respect to the initial system 
(see also iBezanson et al.l (120091 ) for a similar derivation of the scaling relations). We know 
from the shape of the Schechter function for the distribution of stellar masses that a massive 
system (m > M*) accretes most of its mass from lower mass systems and thus the simple 
calculation above makes it very plausible that even though major mergers do occur minor 
mergers are the main driver for the evolution in size and density of massive galaxies. 



3. High resolution simulations of an individual galaxy halo 



We have performed a cosmological N-body/SPH high- resolution re-simulation of an 
in dividual galaxy hal o. The process of setting up the initial conditions is described in detail 



m 



Naab et al.l (120071 ) and is briefly reviewed. W e have re-run galaxy A at 200 3 particles 



resolution using a WMAP-1 JSpergel et all l2003f ) cosmology with a Hubble parameter of 
h = 0.65 (= H =100h kms^Mpc" 1 ) with a 8 =0.86, f b = Q b /Q m =0.2, fi =0.3, and A =0.7. 
To re-simulate the target halo at high resolution we increased the particle number to 1.6 x 10 7 
gas and dark matter particles within a cubic volume at redshift z = 24 containing all particles 
that end up within the virialized region (we assumed a fixed radius of 0.5 Mpc) of the halo 
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at z — 0. The tidal forces from particles outside the high resolution cube were approximated 
by increasingly massive dark matter particles in 5 nested layers of decreasing resolution. The 
galaxy was not contaminated by massive boundary particles within the virial radius. 



The simulation was performed with GADGET-2 (ISpringell 120051 ) on Woodhen at the 
Princeton PICSciE HPC center using a total of 177,000 CPU hours on 64 CPUs. We used 
a fixed comoving softening until z = 9, and thereafter the softening, e.g. for the stars, 
remained fixed at physical 125pc. The mass of an individual stellar particle is 1.3 x 1O 5 M 
and we spawn two stars per SPH particle. 

Star formation and feedback from supernovae was included using the sub-grid multiphase 
model of (jSpringel &: Hernquistll2003l ). We require an over-density contrast of A > 55.7 for 
the onset of star formation to avoid spurious star formation at high redshift. The threshold 
number density for star formation is n t h re sh — 0.205 cm -3 and the star formation time-scale 
is i* = 1.5/j, -1 Gv rs. We also includ ed an uniform UV background radiation field peaking at 
at z ~ 2 — 3 



sec 



Naab et all 120071 ). 



At present the galaxy has a total virial mass of 1.9 x 1O 12 M and a total stellar mass 
of 1.5 x 10 11 M Q . The ratio of central stellar mass to halo mass is about a factor of two 



larger than predicted from gravitational lensing studies (IMandelbaum et al.ll2006l ) , however. 



it is c omparable to some recent pred ictions derived for the Milky Way halo ( e.g. iXue et al. 
20081 . see however iLi &: White! 120081 ) . The central stellar co mponent resembles an elliptical 
galaxy with properties very similar to the results presented in iNaab et al.l (120071 ) and in this 
letter we only focus on particular aspects of the assembly and size evolution of the stellar 
component. 



4. Results 

During the assembly of the central galaxy we have separated the stars within a fiducial 
radius of 30 kpc in fixed physical coordinates into stars that have formed in-situ from gas 
within the galaxy and stars that have formed outside this radius and were accreted later-on. 
The mass assembly history of the in-situ and accreted components of the stellar system are 
shown in Fig. [TJ The early mass evolution at z > 2 is driven by the assembly of in-situ 
stars with a decreasing contribution towards z ~ 0.7. Below this redshift only few stars are 
formed within 30 kpc. The final 20% of stars are added thereafter by accretion of systems 
formed outside the main stellar system at radii larger than 30 kpc. 

The upper panel of Fig. [2] shows the density profiles of the in-situ stars at redshifts 
z=5,3,2,l,0. Between z=5 and z=3 the central galaxy is still building up from gas flows 
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feeding the central region of the galaxy directly, forming a concentrated stellar system. The 
in-situ central stellar densities decrease by more than an order of magnitude towards lower 
redshifts. The spherical half-mass radii of the in-situ ste llar component, show that the 
in-situ system is very compact (see also IJoung et al.l 120091 ) at z = 3 (0.5 — 0.6 kpc) and 
its size increases by about a factor of four (w 2 kpc) until z=0. In the bottom panel of 
Fig. [2] we show the density profiles for the stars that have formed outside 30 kpc and have 
been accreted later-on. This component is more extended at all redshifts and has a shallower 
density profile. Its central density stays almost constant at ~ lO lo M kpc~ 3 while the density 
at larger radii subsequently increases towards z = 0. The half- mass radius of this component 
is significantly larger than for the in-situ stars (> 3 kpc). The central part of the galaxy is 
always dominated by in-situ stars whereas at redshifts below z ~ 2 — 3 and at radii larger 
than « 2 — 3 kpc the system is dominated by accreted stars. 

In Fig. [3] we show the time evolution of the edge-on projected half mass radius of stars 
in the central galaxy within the central physical 30 kpc as a function of time. At z = 3 
the stellar system resembles a compact disk-like or bar-like object with a peak ellipticity of 
e = 0.65 and a size of ~ 0.3 — 0.7 kpc at z=3. Thereafter its size increases by a factor of 
~ 3 — 4 to its present value of 2.7 kpc. We also plot the proje cted half-light radii in the 
rest frame K- and V-band using the stellar population models of iBruzual fc Charlotl (120031 ) 
assuming solar metalicity. In general the half-mass radii trace the half-light radii even at 
larger redshifts reasonably well. 

The K-band rest-frame surface brightness profiles for edge-on projections at z=0 and 
z=3 are sho wn in Fig. H]in combinat ion with the best fitting Sersic profiles, using the fitting 
procedure of iNaab fc Trujillol (120061 ) excluding the central three softening lengths. At high 
redshift the system is very compact, r e = 0.69 kpc, and has a moderate Sersic index of 
n ~ 2.3. This is in agreement with the system being flattened and disk-like. At low redshift 
the system is more extended r e = 2.4 kpc, and its Sersic index has increased to n = 3.3. The 
galaxy is slightly more compact t han typical SPSS early- ty pe galaxies at this mass but lies 
within the observed distribution (IShen et al.ll2003l . see also iFranx et al.l 120081 ) . The errors 
given in the figure are bootstrap errors for a fixed projection. As we have shown before the 
evolution in surface brightness is mainly driven by an evolution in surface density and not 
by stellar evolution. 



At z=3 the system has a total stellar mass of M = 5.5 x 10 
radius r cS = 0.69 kpc and a corresponding effective density of p c s = 0.51 
1.6 x 10 10 M Q . The projected stellar line-of-sight velocity dispersion is cr c g- 



M P) with an effective 

1 K 1\/T I ( A /Q^3 \ _ 



240 kms" 1 . 

The corresponding values at z=0 are M = 15 x 1O 1O M , r cS = 2.4 kpc, p e ^ = 1.3 x 10 9 M Q 
and (TefT ~ 190 kms -1 , which is a typical dispersion for early-type galaxies at this mass 
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( IBender et al.lll992l ). From z=3 to z=0 the system accretes about 5.5 x 1O 1O M (see Fig [I]) 
and we can assume for the above scaling relations rj = 1 and e << 1 (which is a reasonable 
assumption as most mass is accreted in very small systems and we found the most massive 
merger since z=3 with a mass ratio of 8:1). The z=0 values of the simulated galaxy are close 
to the simple prediction, however, the size increase as well as the decrease in density and 
dispersion are more moderate. This is however expected as the real evolution of the system 
is more complex and there is non- negligible in-situ star formation between z=3 and z=l. 
Still, the simple scaling relations for stellar accretion represent the evolution of the system 
from z=3 to z=0 very well. However, observations of more massive ellip ticals than the one 
presented here indicate an even stronger size increase (IFranx et al.l 120081 ). This effect can be 
expected if the assembly of more massive galaxies is even more dominated by minor mergers 
and stellar accretion. 

In particul ar, the drop in veloci t y disp ersion is in qualitative agreement with first direct 
observations by lCenarro fc Trujilld (120091 ). Recently, iBezanson et al.l (120091 ) have reported 
a relatively weak evolution of the density within fixed lkpc of only a factor 2-3. This also 
is in qualitative agreement with our simulation which shows a decrease of only a factor 1.5 
from z=2 to z=0. We also note that the stellar population of the system is already evolved 
at high redshift. At z = 2 the galaxy has a stellar mass of 7.5 x 1O 1O M and a local star 
formation rate of only w 6M Q yr _1 and an average stellar age of 1.6 Gyrs. 



5. Conclusion & Discussion 



In this paper we show that the observed size and density evolution of massive spheroids 
agrees with what is to be expected from a high resolution cosmological simulation of a sys- 
tem which grows at late times predominantly by minor mergers and accretion of stars. We 
can successfully apply simple scaling relations derived from the virial theorem to demon- 
strate that that the size increase and decrease in density and velocity dispersion is a nat- 
ural consequence of mas s assembly by much less massive stellar systems and accretion 
( jCenarro fc Trujilld 120091 ) and cannot be explained by mass assembly histories dominated 
by major stellar mergers. 

In the simulation, a fi rst phase (6 > z > 2) dominated by in-situ star formation from 
inflowing cold gas (see e.g. iDekel et al.ll2009l and references therein) produces a massive and 
dense stellar system with sizes r\/2 < 1 kpc. This phase of the formation of the cores of 
ellipticals is followed by an extended phase (3 > z > 0) with little in-situ star formation 
but significant accretion of stellar material. This material can be stripped at larger radii 
and increases the size of the system with time. At the same time the central concentration 
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is reduced by dynamical friction from the surviving cores (see lEl-Zant et al.ll200ll ). The 
apparent size increase is caused by the initial dominance of the in-situ component being 
heated and, at larger radii, overshadowed ultimately by the accreted stars. From z=3 to z=0 
the effective radius of the system increases by a factor of 3.5 with a decrease in the effective 
density of more than an order of magnitude and a decrease in velocity dispersion of 20%, 
in good agreement with predictions from simple scaling relations for the accretion of minor 
mergers. 



Genel et al.ll2008bl: 


Fakhouri & Ma 


g- 


Bundv et al. 


2009) on the 



mass assembly mechanisms of early-type galaxies and their dark matter halos demonstrate 
the significance of minor mergers. In addition, due to the shape of the mass function, massive 
systems at high redshift are more likely to experience minor mergers, than lower mass galax- 
ies. If the size evolution is in general driven by minor mergers we would expect a differential 
size increase, e.g. more massive high redshift systems grow larger than lower mass systems. 
At the same time minor mergers do also play an important r ole for the gravitational heating 



of halo gas, thereby suppressing the formation of new stars (1 Johansson et al.ll2009bl ). 



A picture of a two phase formation process for massive spheroidal galaxies has a number 
of virtues. In the first dissipative phase at high redshift stars form quickly and build the 
compact progenitor of present-day ellipticals. In fact, it seems of minor importance if the gas 
is fu nneled to the center through streams or merg e rs of extended ga s dominated disks (see 



e.g. iRobertson et al.ll2006l ; I Johansson et al.ll2009al ; iDekel et al.l 120091 ) as long as it happens 
on a short timescale. The stars formed at this early pha se are expected to b e significantly 
enriched in alpha-elements as expected from observations (IThomas et al.ll2005l ) and form the 
compact core of the elliptical galaxy. Later on metal-poor stars from smaller systems are 
accreted and form the halo of the galaxy resulting in the observed metalicit y gradient. This 



i nside out formation scenario is also made plausible by recent observations of iBezanson et al. 



( 120091 ). There is an important test of the picture presented in this paper. If correct, then the 
outer parts of massive giant ellipticals will tend to be old, blue, metal-poor and relatively 
uniform from galaxy to galaxy since they are all composed essentially of the debris from 
tidally destroyed accreted small systems. 



The simulation was performed at the Princeton PICSciE HPC center. This research 
was supported by the DFG cluster of excellence 'Origin and Structure of the Universe'. We 
thank Marijn Franx, Pieter van Dokkum and Ignacio Trujillo for helpful comments on the 
manuscript. 
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Fig. 1. — Mass assembly history of the stellar system (squares) separated into stars made 
in-situ (open diamonds) in the galaxy and stars formed outside the galaxy that have been 
accreted (stars) later on. At high redshift (z > 2) the system assembles by the formation of 
in-situ stars, at low redshift [z < 1) accretion is more dominant. 
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Fig. 2. — Density profile of the stars formed in-situ in the galaxy (upper plot) and of stars 
formed outside the galaxy and then accreted later-on (bottom plot) at redshifts z=5,3,2,l,0 
(red, orange,blue, green, black). The spherical half mass radii r\n are indicated by the dashed 
vertical lines. The shaded area indicates the gravitational softening length. 
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Fig. 3. — Time evolution of the edge-on circular projected stellar half-mass radius within 
fixed physical 30 kpc (black diamonds). From z = 3 — 2 to z = the size increases by 
a factor of ~ 3 — 4. We also show the evolution of the rest-frame K-band (red cross) and 
V-band (green triangle) half-light radius. Outliers indicate minor merger events, e.g. the 
most massive (8:1) merger since z=3 at z=0.3. 
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Fig. 4. — K-band surface brightness profiles assuming solar metalicity at z = 3 (black dots) 
and z = (red dots). The black and red lines indicate the best fitting Sersic profile. At high 
redshift the galaxy has a higher central surface brightness and is more compact. 



